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Past work has shown that ions can pass through a membrane more readily in one 
direction than the other. We demonstrate here in a model and an experiment 
that for a mixture of small and large particles such asymmetric diffusion can arise 
solely from an asymmetry in the geometry of the pores of the membrane. Our 
deterministic simulation considers a two-dimensional gas of elastic disks of two 
sizes diffusing through a membrane, and our laboratory experiment examines the 
diffusion of glass beads of two sizes through a metal membrane. In both experiment 
and simulation, the membrane is permeable only to the smaller particles, and the 
asymmetric pores lead to an asymmetry in the diffusion rates of these particles. 
The presence of even a small percentage of large particles can clog a membrane, 
preventing passage of the small particles in one direction while permitting free 
flow of the small particles in the other direction. The purely geometric kinetic 
constraints may play a role in common biological contexts such as membrane ion 
channels. 
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Asymmetric diffusion appears in studies of ion transport across membranes (0, 0, 0, 0, 0, 0, 
, [l~2l ) , and in the context of osmosis \l4 ). Ion pores in biological membranes 
often act as effective rectifiers, blocking the passage of a particular ion in one direction while 
allowing free flow in the other. Explanations of these phenomena emphasize electrostatic effects 
or conformational changes. For example, a blocking ion might pass part way through a membrane 
pore, and then become bound, preventing passage of subsequent ions. If the pore is asymmetric, 
the membrane can act as a rectifier. Here we demonstrate that asymmetric diffusion can occur 
with no ionic binding, in a fixed pore geometry. Blocker particles can be held in place by purely 
geometrical and kinetic constraints, even though particles are in rapid random motion. 

We present both a laboratory experiment and a computer simulation to demonstrate this 
effect. In the experiment, macroscopic (1 mm) glass beads diffuse through a metal asymmetric 
membrane and random motions are generated by rapidly vibrating the container. In the simula- 
tion, a two-dimensional gas of elastic hard disks of different sizes diffuses through an asymmetric 
membrane. In both cases, nearly complete rectification is easily demonstrated. 
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Results 

Simulations. Figure [I] shows typical snapshots from the hard disk gas simulation. In both 
panels the gas is initially restricted to the right chamber, and then tends toward equilibrium 
as the smaller disks move from right to left through the membrane. But in the lower panel, 
the larger disks soon clog the pores and greatly reduce the rate of diffusion of the smaller disks 
through the membrane. The presence of the smaller disks bombarding the larger disks within 
the pore makes it highly unlikely that the larger disks will be able to emerge from the pore. 

Figure [2] displays the concentration in the left chamber as a function of time for both orien- 
tations of the membrane. The flow entering from the smaller end of the pores (Fig. [T]A) behaves 
as if there were a well-defined fixed diffusion constant, but the flow through the membrane that 
enters the larger end of the pores (Fig. [T]B) drops to near zero, as the membrane seals itself. The 
clogging behavior of this direction of flow is maintained even at lower concentrations of large 
disks. The number of disks in the plateau region may vary from run to run, depending on how 
many small disks are able to cross before the membrane seals. 

The simulated model is a purely energy-conserving Hamiltonian system, yet it can relax to 
a long-lived macroscopic state of disequilibrium, with two subsystems at different temperatures. 
The more rapidly moving disks preferentially pass through the membrane, resulting in a higher 
temperature on the left side. Such a temperature difference is typical of an effusion process (1 15l). 
There is a very slow equalization of temperature through collisions at the small ends of the 
pores, and of concentration through the occasional unblocking of a pore due to fluctuations, but 
the disequilibrium can persist for very long times. 
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Figure 1: Snapshots of a binary mixture of hard sphere disks, initially all on the right side, that 
diffuse through an asymmetric membrane whose pores are smaller on the right end in A and on 
the left end in B. Flow through the membrane is highly reduced with the membrane geometry 
of the lower panel. The larger particles in the gas mixture are slightly too large to pass through 
the pores, and the smaller particles are just small enough to pass. The left chamber is initially 
empty, and the elapsed time is the same for A and B. 
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Figure 2: The number of particles in the (initially empty) left-hand chamber (cf. Fig. [T]) in 
simulation and experiment. In each graph, the upper curves correspond to pores with their small 
end on the right, and the lower curves correspond to the reverse orientation of the pores. The 
upper curves show a large diffusion rate (slope of the curves), while the lower curves illustrate 
that the pores rapidly clog and the diffusion rate approaches zero. The inset in the lower graph 
shows the clogging for three different experimental trials. The right-hand chamber initially had 
1350 disks in the simulation, and 2000 beads in the experiment, which used a membrane of 
thickness 9.53 mm; the data shown are an average of three experiments. In both cases 20% of 
the particles were large and the remainder small. 
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Experiments. The experiment has two species of glass beads that correspond to the two species 
of disks in the simulation, and a brass membrane with asymmetric pores, as illustrated in Fig.[3j 
The results from experiment are qualitatively the same as those from simulation (Fig. [2]). Once 
again, there is free flow through the membrane from the side with the small pore ends, and 
clogging for flow in the other direction. The inset shows that the number of small beads passing 
through the membrane before it clogs can again vary from run to run. The clogging is not 
permanent. Runs longer than 30 min show a leak rate of about 0.1 bead per minute per pore, 
for this number of beads. 




Figure 3: (A) Diagram of the experimental system. (B) Cross-section of a circular pore that 
has a reduced diameter on the right end, too small for the larger particles in the binary mixture 
to pass. (C) Image for 50 beads in a chamber. (D) Image for 250 beads in a chamber. 

Figure H] shows the effect of membrane thickness on the blocking effect. For the thinnest 
membrane, there is no blocking effect, and the flow is even enhanced in the funneled direction. 
This is in accord with the single-species small bead experiment to be mentioned below. The 
membrane must be of a minimum thickness to produce a significant backing up of the blocking 
beads. If there is a large concentration gradient across the pore, the blocking beads will be held 
in place by this gradient. The backing up requirement is central to the effect. 

When only small beads are present in the experiment, there is a bias in the dynamics, as 
Fig. E] illustrates. The small beads pass more readily though the membrane when they enter the 
large end of the pores, and the system approaches a dynamic equilibrium with more beads on the 
side with the small pore ends. This asymmetry can be explained by the dissipative nature of the 
particle- wall collisions, which reduces the normal component of the colliding particle; therefore 
the angle between the particle velocity vector and the sidewall becomes smaller (|16l ). In this 
way the collisions inside the large end of the pore focus a particle trajectory towards the small 
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end of the pore and therefore increase the effective "cross section" . 
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Figure 4: Number of small beads measured to cross the membrane after three minutes, as a 
function of membrane thickness, for the two flow directions (cf. Fig. [3]). The loaded side of the 
membrane initially had a mixture of 1600 small beads and 400 large beads. 



Leaky membrane simulation results. Some insight into the dynamics of the hard disk 
system can be garnered by considering the response to an imposed variation in the number of 
small particles in the chamber on the membrane side with the small pores (left side in Fig. [T]B). 
The resultant macroscopic state of the chamber on the right side, which we will call the "inside" , 
shows a path dependence, as illustrated in Fig. [6l Initially the system had 200 large particles 
in the inside (right) chamber and no large particles on the outside (left) chamber; initially no 
small particles were in either chamber. In Fig. (6JA, the number of small disks in the outside 
chamber is gradually increased at a rate of two disks per time unit, up to 1200 disks. Each disk 
is added with unit speed, at a random position and moving in a random direction. The number 
of small disks on the inside chamber can be seen to follow this increase, as one would expect 
from a simple diffusion process. Then the number of small disks on the outside is reduced at 
the same slow rate, and again the disks on the inside roughly follow. When there is no strong 
concentration gradient across the membrane, the pores conduct freely in both directions. 

In Fig. [6]B, the number of small disks in the outside chamber is again increased at a rate 
of two disks per time unit up to 1200 disks, and then this time the number of disks on that 
side is suddenly reduced to zero; any small disks leaking through from the inside to outside are 
removed from the system. In this case, the instantaneous formation of a large concentration 
gradient across the membrane causes it to clog, and diffusion from the inside initially drops to 
nearly zero, even though the final concentration outside is the same as in Fig. EH- 
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Figure 5: Number of beads measured to pass through a membrane (9.53 mm thick), as a 
function of time, when only small beads were present. Flow in the direction from the large to 
the small end of the pores is enhanced relative to that for the reverse orientation of the pores. 
There was a total of 200 beads; those residing in the pores were not observed. 



The diffusion outward does not, however, remain at zero, because the initial concentration of 
blocker disks inside is small enough so that there is a significant leakage current. Each small step 
in the descending portion of the curve in the lower panel corresponds to a single pore opening 
for a time; a large blocker disk is dislodged by a fluctuation, and a "patch clamp" observation 
of that pore would show a current of small disks flowing outward until the pore is resealed by 
another large disk entering it. Eventually, the concentration of small disks on the right becomes 
too small to maintain the clogged membrane, and the concentration curve on the right returns 
to a simple diffusive decrease. Because there are only nine pores in this simulation, the length 
of the plateau corresponding to the clogging is variable, but the general shape of the knee is 
rep eat able. 

The possible biological relevance of this clogging effect to membrane rectification is evident. 
Possibly the simple rectification described here could have played a role in tide-pool scenarios 
of the origin of life. In our simulation, repeated sawtooth driving on the outside will maintain a 
concentration gradient indefinitely, even with a leaky membrane. Similarly, one can imagine that 
repeated evaporation and flooding of a primeval tide-pool may have provided the conditions for a 
sustained gradient across the membrane of an early vesicle, enough to drive a proto-metabolism 
0). 

Discussion 

If flows across a membrane are governed by the ordinary diffusion equation, they will reverse 
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Figure 6: Time-dependent variation of concentration of small disks (solid lines) in the chamber 
on the right (the "inside" ) in Fig. [T]B, in response to variation of number of small disks in the 
left-hand ("outside") chamber (dotted lines). The number of large disks in the inside chamber 
is 200 and is zero in the outside chamber. A When the number of disks outside is linearly 
increased to 1200 and then decreased linearly at the same rate, the number of disks on the 
inside adiabatically follows. B When the number of small disks outside is linearly increased as in 
A and then suddenly decreased to zero, the membrane pores clog as a consequence of the sudden 
strong cross- membrane gradient, and the number of small particles on the inside remains nearly 
constant for a long time. 



exactly, if the concentration gradient across the membrane is reversed, regardless of the geometry 
of the pores. Implicit in the use of the linear diffusion equation is the point-particle assumption. 
Particles do not interact, and solutions to the equations of motion are the same up to a factor 
for a dense gas or a dilute gas, both in the bulk and in the pore. 

If the diffusing particles have a size comparable to the pores through which they are mov- 
ing, the linear diffusion equation is no longer applicable. In the dilute limit, flows will still be 
symmetric under reversal of the concentration gradient, for any pore shape. But when concen- 
trations are raised and particles interact with each other and the walls of an asymmetric pore, 
the symmetry is broken, and the diffusion rates from one side or the other can differ. 

Flux ratio theorems indicating a symmetry in diffusion across membranes have been proven 
for a variety of systems, e.g. with linear intra-membrane trapping kinetics for the diffusing par- 
ticles (3, Jjjl) and even for multicomponent systems such as the one studied here (0). However, 
these theorems fail here because the trapping kinetics do not include the kinetic constraints felt 
by the spatially extended particles interacting with pores. 

At first examination asymmetric diffusion might seem counterintuitive, at odds with the 
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Second Law. But there is no contradiction, as this is a nonequilibrium effect, being driven 
here by increased concentration on one side. At equilibrium, detailed balance will hold, and we 
expect fluctuations around equilibrium to be symmetric. But a clogged membrane is held away 
from equilibrium, in a long-lived metastable state. A concentration gradient across the pore is 
required to keep blocking particles in place, but this gradient can be made arbitrarily small by 
lengthening the pore. 

Asymmetric diffusion could be relevant to biological situations; rectification can be generated 
in isothermal conditions, driven by only a small concentration gradient. Asymmetry has long 
been noted in movements of ions through channels in membranes, where the asymmetry is often 
described as a form of rectification |l|). Although ion channels can be complicated devices, with 
voltage sensors and internal conformational changes, a particularly simple pore design is the 



inward-rectifier potassium channel; for a review see Lu (|2ll ). Early modeling of this channel 
described the outward flow of K + ions through the channel as being blocked by lar ger ions, 
which are present on the inside of the cell and cannot pass completely through 



Recent detailed imaging (|26|) of a potassium channel shows an asymmetic pore, narrow towards 
the outside and wider towards the inside, roughly like our model pores represented in Figs. [3]B 
and[Q 

The early rectification models of Hille and Schwartz (jiil ) explicitly included the requirement 
that a blocker back out of the channel to re-enable flow. However, these and other models 
emphasize electrical effects, in particular intra-pore binding sites, to keep the blockers in place 
in the pore. Electrostatic interactions among ions and possible binding sites will of course be 
very important, but we have demonstrated that an asymmetric geometry alone can also lead to 
rectification. We suggest that the clogging effect described in this paper could play a role in 
maintaining blocking and rectification. The inclusion of this effect may help clarify discussions 



of "multi-ion channels" (|23|, 1271 ). Geometry- induced asymmetric diffusion could occur in any 
physical context that includes the interaction of shapes in restrictive geometries, as exemplified 
by our experiment and its membrane model. It is natural for large objects to get stuck and to 
block passage of following objects; all that is required is a certain density and a flow direction, 
no matter how small the flow. The simplicity of the conditions generating the effect suggest 
that it might occur quite commonly in biological systems, where membranes are ubiquitous and 
molecules of all shapes and sizes are closely packed (j28l ) . Large changes in forward and backward 
diffusion rates can be effected by small changes in pore geometry, or object shape. The simple 
clogging phenomenon arising from asymmetric pores suggests a possible evolutionary path from 
a purely mechanical rectification that might occur in a pre-biotic vesicle, to the sophisticated 
membrane ion channels present in living systems today. 

Materials and Methods 

Experiment. A parallelepiped container is made of standard microscope slides for the top 
and sides, and cut slides for the ends. The bottom of the container is brass with the interior 
surface machined with a pattern of 2 mm wide divots to deflect beads from vertical motion, to 
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help randomize the motion. The container is divided into two chambers by a brass membrane 
whose thickness D was varied (Fig.[3H)- Each membrane has 16 circular pores with a long large- 
diameter section (whose length depends on the membrane thickness) and a concentric short 
small-diameter constriction (whose length is fixed, 1 mm), as shown in Fig. [3]B. 

Particles of two different sizes are obtained by repeated sieving of soda lime glass beads 
(sphericity > 95%) with a broad size distribution (Ceroglass GSR). The "large" beads (average 
mass 7.6 mg) have diameters 1.70-1.85 mm, small enough to fit into the large-diameter end 
of a membrane pore but too large to pass through the narrow end of a pore. The "small" 
beads (average mass 4.2 mg) have diameters 1.40-1.55 mm, small enough so they can easily pass 
completely through the pores, but large enough so that they have to stay in single file in the 
large-diameter portion of the pore (see Fig. [3]B). 

In each run the chamber on one side of the membrane is loaded with beads, and then 
the container is oscillated vertically at 30 Hz with an acceleration amplitude 9.5c? (g is the 
gravitational acceleration). To the eye, the bead motion is vigorous and random. Pictures such 
as those in Fig. [3]C and D are taken at one-second intervals (1/2000 s exposure time) with a 
digital imaging system (Kodak Motion Corder Analyzer SR 1000) synchronized with the driving 
frequency. 

The number of small beads passing through the membrane into an initially empty chamber 
is counted as a function of time using an image processing algorithm that is accurate to within 
one or two beads for up to about 80 beads, as in Fig. [3]C. For a larger number of beads, as in 
Fig. [3p, the algorithm undercounts the number of beads due to eclipsing and grouping effects, 
but a correction factor obtained by weighing the beads and checked by direct hand counting 
yields counts accurate to 5%. 

Hard disk simulation. The simulation uses a two-dimensional container divided by a mem- 
brane into two chambers, each 3 units high by 2 units wide, as shown in Fig. [TJ The membrane 
has nine tapered pores, 0.04 units wide on the large end and 0.02 units wide on the small end. 
There are two species of elastic hard disks, "small" disks (diameter 0.018) and "large" disks 
(diameter 0.022). As in the experiment, the small disks can pass through the pores in either 
direction, while the large disks cannot pass through the narrow end of a pore. There are initially 
1350 disks on the right side, 20% of which are the large disks; the left side is initially empty. 

All disks have the same mass and have initial unit velocities uniformly distributed in random 
directions. A period of at least 10 time units is allowed to pass before the pores are opened, to 
allow the velocity distributions to thermalize; each simulation is subsequently run for 500 time 
units. The work of Forster et al. (0) shows that the dynamical entropy of such systems is high, 
so that the Maxwellian distribution is rapidly achieved, and preserved even in very restrictive 
geometries. 
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